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Abstract 

By means of simple models in a fiat spacetime manifold we examine some of the issues that 
arise when quantizing interacting quantum fields in multi-metric backgrounds. In particular we 
investigate the maintenance of a causal structure in the models. In this context we introduce 
and explain the relevance of an interpolating metric that is a superposition of the individual 
metrics in the models. We study the renormalisation of a model with quartic interactions and 
elucidate the structure of the renormalisation group and its implications for Lorentz symmetry 
breakdown. 
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1 Introduction 



In standard General Relativity the structure of the space-time manifold is expressed in terms of 
a metric field. This metric provides a lightcone structure that determines causal relationships 
between events in the manifold. Modified gravity theories involving more than one metric have 
been proposed as potential explanations of dark matter and dark energy effects in cosmology 
[1, 2, 3, 4]. They are also relevant to models of massive gravity[5, 6]. A comprehensive report 
on work in this area can be found in reference [7]. While these approaches are for the moment 
conjectural they are of considerable interest in themselves. 

An important feature of such theories is that they may lead to a violation of Lorentz invari- 
ance. There have been many careful studies of the breakdown of Lorentz invariance [8, 9, 10] 
with the conclusion that observations constrain it to be a very small effect indeed [11]. Never- 
theless it is interesting to consider what issues of principle are involved in dealing with quantum 
field theory in the context of a multi-metric space-time without any a priori constraint on the 
relationship of the various metrics except for the general requirement of maintaining a causal 
structure for the theory. These issues of causality and multiple lightcone structure appear in in 
other contexts [12], [13, 14, 15], [16]. 

In this paper we are concerned not with the gravitational dynamics implied by such theories 
but with the behaviour of quantum fields in the presence of more than one space-time metric. We 
therefore simplify the problem and consider a flat space-time with several metrics each of which 
is associated with a particular quantum field. We investigate the constraints resulting from a 
requirement of causality by studying some simple scalar field theory models. In the simplest 
case we have a model with two fields and two metrics. We find that if we wish to enforce 
causality when there are interactions between the fields, it is important that the two metrics 
support lightcones that overlap and share spacetime vectors that are timelike in both metrics. 
A useful and significant construct that plays a role in elucidating the relationship between the 
two metrics is that of the interpolating metric which is built from a linear superposition with 
positive weights of the two metrics. As the weights of the superposition vary the interpolating 
metric passes from one metric to the other and may at an intermediate stage become singular. 
The avoidance of this singularity is important for the dynamics of the theory and leads to the 
restriction indicated above requiring the two lightcones to intersect by sharing a set of spacetime 
vectors that are timelike in both metrics. We indicate briefly generalisations to cases with more 
than two metrics. 

For a model in which the scalar fields experience quartic interactions we investigate the 
renormalisation of the perturbation series. The interpolating metric plays a crucial role in the 
evaluation of the of the Feynman diagrams of the theory reinforcing the restrictions on the 
lightcones required for a causal structure. The individual metrics are renormalised as part 
of the calculation and are involved through the properties of the interpolating metric in the 
renormalisation group flows describing the evolution of coupling constants and masses in the 
theory. 



The fundamental assumption we make is that there exists a flat space-time manifold the events 
in which are specified by a set of coordinates {^i — 0,1,2,3} that are arbitrary under 
linear transformations. On the manifold are a set of N metrics gfy, {i = 1, 2, • • • , N} that are 
independent of . Each metric specifies a lightcone on the manifold through the equation 



That eq(l) does determine a lightcone is ensured by constructing g^l from an associated vierbein 
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Figure 1: Examples of pairs of lightcones (in 2+1 dimensions) that share spacelike surfaces. The 
first two examples represent situations in which there are spacetime vectors that are timelike in 
both cones. In the third case each lightcone is seen as spacelike from the other. 

Of course the vierbein is arbitrary under Lorentz transformations, L a b , 

e W\^L\eW\ . (3) 

With each vierbein we associate a set of coordinates y^ a = e (<) V . We also have the inverse 
relation . The space-time volume elements are related by 

d 4 y (i) = nWfa , (4) 

where 

fiW =dcte« a M . (5) 
A plane wave in space-time has the form 

f q (x) = e-^ . (6) 
Viewed from the vierbein frame (i) the four vector is 

p W« = e W^ . ( 7 ) 

Obviously 

q^=P {t) a y {l)a , (8) 

and 

9 {i)ftv q^ = VabP {i)a p {i)b ■ (9) 

Our aim is to investigate models of quantum field theory in which each lightcone is associated 
with a set of fields that transform in a conventional fashion under a Lorentz group associated 
with that lightcone but which may also interact with fields associated with other lightcones. 
This is a scalar field version of a model investigated by Colemam and Glashow [17] in relation 
to neutrino mixing. There are then two issues. The first is that there is in general no overall 
Lorentz invariance for such a theory. There may however be circumstances, as suggested by 
Coleman and Glashow [8], in which a reduced Lorentz invariance survives. 

The second issue is that we require restrictions on the relationship between the lightcones 
in order that there can be a causal evolution of the fields. A natural way of achieving this 
is to require that there exists a foliation of the (flat) space-time manifold that is space-like 
with respect to each of the lightcones simultaneously. An example of the possible relationships 
between two lightcones that respect this constraint are illustrared in Fig 1. 

We will identify the coordinate that labels the leaves of the foliation, as the time variable. 
Of course if there is one such foliation we expect there will normally be an infinite set. This 
will become clearer in the example of two lightcones described in section 3. When we come to 
consider the dynamics of interacting quantum field theory models in more detail we will be led 
to impose further restrictions on the relationship between the lightcones. 

In order to emphasize the importance of the space-time manifold we will formulate the theory 
so that it is invariant under general linear transformations of the coordinates on the manifold. 
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3 Bi-Metric Geometry 

The basic ideas can be understood by considering a model in which there are two metrics. We 
simplify the notation and introduce metrics and g M „ that are constructed from two vierbcins 
e afj, and e af i . 

g^u = Vabe a M e b v . (10) 
This guarantees that the two surfaces 

g^vx^x , 

g^afx" = , (11) 

are indeed lightconcs of conventional form. The relationship between the two lightcones deter- 
mines the nature of the breakdown of Lorentz invariance in the theory. We have 

ft = dete Q M , 

f2 = dete Q M . (12) 
The inverse vierbeins are e M a and e M a and satisfy 

p a A a 

b — °b ! 

e\e\ = 5t . (13) 
We can construct coordinates appropriate to each vierbein, 

y a = eV" . 

V a = e°X . (14) 

= , 

= e\Tf ■ (15) 
a ^L\y b . (16) 



The inverse transformations are 



Under a Lorentz transformation L a b 



Now 

e a ^L\e\ = e a v K% , (17) 

where therefore A" is also a representation of the Lorentz group that preserves the metric g^ u . 
Under the same Lorentz transformation it follows that 

-> A" v x v . (18) 

Similarly for the vierbein e afJi , we have associated coordinates y a where 



and under a Lorentz transformation L a b 



with 
where 



(19) 

L\v b , (20) 
K» v x v , (21) 



e%^L\e\ = e a v K\ . (22) 

It follows that K u is a representation of the Lorentz group that preserves the metric g^ v . 
It is obvious that in general v and v arc not the same representations of the Lorentz 
group since they preserve different metrics on space-time. As indicated in section 2 there can 
be circumstances in which the two versions of the Lorentz group share a subgroup [8, 17]. This 
subgroup will then provide a reduced Lorentz invariance of the overall theory. 
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3.1 Foliation and Linear Transformations 



The coordinate space of the flat manifold is identical with the tangent space (at the origin). In 
order to set up a foliation of the coordinate vector space we choose a cotangent vector n M and 
define a subspacc X by the restriction 

= . (23) 
We then introduce a vector m M not lying in X which satisfies 

n^m" = 1 . (24) 
A general vector can be expressed in the form 

x» = rm^ + j/ 1 , (25) 

where y M lies in X . It follows that 

n^ = r . (26) 

The subsets of coordinate space each labelled by a value of r, constitute a foliation. 

It is useful to parametrize linear transformations on coordinate space in a manner consistent 
with the foliation. First we introduce general linear transformation M 11 v on the subspace X 
which satisfies 

n^M\=n v , (27) 

and 

M\y v = y» . (28) 

This transformation maps X into itself and leaves r unchanged. We follow with a shear trans- 
formation S(yY v that has the form 

S"„=J£-v"n„ , (29) 

where is a vector lying in X, that is n^v^ = . It follows that 

Si* v (v)x v = x" - tv» . (30) 

Clearly 5 M v (v) leaves t unaltered but shifts the origin of coordinates lying in the subspace X 
by an amount proportional to r . We have also 

{S-\v)Y v = S\{-v) . (31) 

If we assume coordinate space has n dimensions then X has n—1 dimensions and the linear 
transformation M M „ has (n — l) 2 parameters. The shear transformation has n — 1 parameters 
so the combined transformation (S(v)M) 11 v has n 2 — n parameters. If we now apply an overall 
scale transformation, uj5^ this introduces one further parameter, bringing the total to n 2 — n+1 . 
This procedure is the most general transformation that preserves the the foliation determined by 
n M and . The remaining n — 1 parameters needed to complete the parametrization of the a 
general transformation may be chosen to be the velocities of the n—1 appropriately constructed 
Lorentz boosts along n—1 directions within the the foliation subspacc X. 

We emphasise that although the resulting parametrization of the linear transformation on 
coordinate space is perfectly general for transformations in a neighbourhood of the identity it 
will not cover the whole general linear group. It is however the appropriate parametrization for 
our purposes. 
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3.2 Relationship between Metrics 

The relationship between two metrics can be associated with a general linear transformation of 
the coordinate vector 

a:" -> at* = A» v x v , (32) 

and 

9^ -> 9^v = A a ^A T u g aT . (33) 

The construction of metrics from vierbeins is an example of this procedure. 

We can now see that it is convenient to parametrize v in the following way, 

A\=uA\S\M\ , (34) 

where A M v is a general Lorentz boost in a representation that leaves g^ invariant. This contains 
the remaining n — 1 parameters needed to complete the total of n 2 for a n dimensional general 
linear transformation, that is connected to the identity. The advantage of this completion is 
that it leaves the relationship between the two metrics unchanged. If we have identified X as the 
(n — 1) dimensional subspace of displacements that are spacelike with respect to both metrics 
then the final Lorentz boosts may be omitted without losing generality. 

Finally we note as indicated above, that the general vierbein (or n-bein) can be constructed 
by the same sequence of steps, the starting metric being the standard Lorentz metric in this 
case. 



3.3 Interpolating Metrics 

A further construct that is important in reflecting aspects of the relationship between two 
metrics is that of the interpolating metric. It turns out to play a crucial role in elucidating the 
causal structure of interacting quantum field theory with two or more lightcones. We define two 
versions. The first is g^(u) where < u < 1 and 

gr{u) = u<r + {i.-u)r v , (35) 

with its inverse g^iu) . The second is g^u) where 

9nA u ) = (1 - u)g^ v + ug^ v , (36) 

with its inverse g^ v (u) . These two interpolating metrics are not inverses of one another but do 
satisfy the relation 

g^(u)=g^~g aT (u)g TU • (37) 

It follows that 

det g^ (u) = det g^ det g ar (u) det g TV , (38) 

and therefore that if one of these interpolating metrics becomes singular so does the other. 
It turns out that g^"(u) is more directly related to the dynamics of the field theory but the 
behaviour of guriu) is more easily interpreted in terms of the the relationship between the 
coordinate lightcones. 

The outcome of a detailed investigation is that if the two lightcones overlap so that there are 
coordinate vectors that timelike with respect to both lightcones then the interpolating metrics 
remain nonsingular for < u < 1 . When no such vectors exist then the interpolating metrics do 
become singular for u in this range. The transitional situation arises when the lightcones of the 
original metrics touch one another along a shared light ray. We illustrate these circumstances 
with a simple example in section 5. A more detailed analysis is presented in appendix A. 
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4 Simple Quantum Field Theory Model 

We construct a simple quantum field theory involving a real scalar field 4>{x) that we associate 
with the vierbein e a . It has a Lagrangian density 

C 4> {x)= l -{g^d l ,cj>{x)d v <P{x)-m 2 ^{x)) . (39) 
The action is given by 

= [ d^xttC^x) . (40) 



Clearly at this stage the action is invariant under both Lorentz transformations of the vierbein 
e a and general linear coordinate transformations. In the notation of subsection 3.1 we choose 
X to be a spacelike subspace of the metric and set — (1, 0, 0, 0, ) and m M = (1, 0, 0,0,) . 
The coordinate x° is the time variable labelling the leaves of the foliation and the coordinates 
x = (a; 1 ,^ 2 ,^ 3 ) are coordinates on the spacelike leaves. 

The quantization is achieved by introducing the field, n(x) 7 conjugate to the field (j)(x) in 
the standard way 

<x) = = Slg^d^x) - n (g m d Hx) + 9°%^)) . (41) 

We complete the procedure by imposing the equal time commutation relations 

[<f>(x°,x),n(x°,x')] = i«J(x-x') , 

[^x),^ ^')] = , 

[7r(x ,x),7r(z°,x')] = . (42) 
The Hamiltonian generating time development is 

= J d 3 x{ir(x , X ){do<P(x , X ))-n£4x ,x)} . (43) 

That is 

= J d 3 x|^^ (tt(x) - n(g 0l d^))) 2 - l^di^dMx) + ^m 2 O0 2 (x)j . (44) 

It is straightforward to check that this Hamiltonian yields a time development for the system 
that is the standard equation of motion derived directly from the Euler-Lagrange equations, 
namely 

g^d^d v 4>{x)+m 2 tf{x) =0 . (45) 
We can obtain solutions for <p(x) as superpositions of plane-waves 

=/,(*)= e-*"*" , (46) 

where q satisfies the mass-shell condition 

gTq„qu = m 2 . (47) 

It is implicit that we select the branch of the mass-shell condition corresponding to positive 
energy in the appropriate vierbein frame. Note that because of our definitions the velocity of a 
particle wave-packet is v = —dqo/dq . 

If f(x) and f(x) both satisfy eq(45) then the current 

J?(x) = XlgT {f*{x)d v f(x) - f(x)d v f'*(x)) , (48) 
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satisfies 

d^(x) = . (49) 

This permits us to define the time- independent scalar product of two wave functions f(x) and 
f'(x) to be 

(/', /) =Mj d 3 x {/'* (x). 9 °^/(x) - f(x)g^dj / *(x)} , (50) 

We have then 

(f g >,f q ) = 2Slg 0v q v (2w) 3 8 ( < l -< l ') , (51) 

where q is the spatial part of . 

As a superposition of plane waves the field <f>(x) can be put in the form 

= / 2^CV) 3 (/g(X)a(g) + f;{xWiq)) ' (52) 
The mode operators a(q) are given by 

a{q) = {f q ,4>) 

= iSlJ d*K{r q {x)g 0v d v 4>{x)-<t>{x)g 0v d„r q (x)) 

= ijd 3 X f*(x)(n(x)-mg »q„cj>(x)). (53) 

Similarly we have 

a\q) = -i j ' d 3 xf q {x){Tr{x)+ing 0v q v <j>(x)) . (54) 

The commutation relations in eq(42) imply 

[a(q),a(q')} = , 
[aHq),a^q')] = , 
[a(q),aHq / )] = 2ft. 9 %,(27r) 3 <5(q - q') . (55) 

The Hamiltonian can be expressed in terms of the mode operators (after normal ordering) 
in the form 



d 3 q 
r)3f>2 3 °>V 



H <t> = J ( 27r) :u,7 WiH'lhi('l) ■ ( r >(>) 
The Feynman propagator is 



G F (x-x') = -i(0\T((/){x)<f){x'))\0) . (57) 

It satisfies the equation 

Wd li d v + m 2 )G F (x-af) = ~5(x-J) , (58) 
and hence can be put in the form 

1 f d 4 a e -iq^ 

Gp(x — x 1 ) — / 77r~TT 5 • (59) 

tK ' ft J {2n) i g^q lx q v -m 2 +ie V ; 

This analysis is no more than standard elementary quantum field theory (QFT-101!) written 
out for general (linear) coordinates. 

We now introduce a second real field ip(x) associated with the vierbein e afl together with a 
Lagrangian density 

■ (60) 
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The associated action is 

S$ = J d i x£^,(x) . (61) 

We can carry out an analysis entirely parallel to that for the field <j)(x) . The field conjugate to 
ip(x) is uj{x) where 

- s(MW) = ^"WW = " (9 m dom + g 0i dMx)) . (62) 

The standard equal time commutation relations are imposed on ip(x) and uj(x) . The equation 
of motion for ip(x) is 

rW(a;) + mVW = , (63) 
and the mode decomposition for tjj(x) in the form 

= 1 2n/ q %)M x »® + J « (x)hHq)) ■ (64) 

The wavefunctions f q (x) are 

f q (x) = e-^ , (65) 
where now q satisfies the mass-shell condition 

g» v qii q v = m 2 . (66) 

The appropriate definition of the scalar product between wave functions has the form 

(/',/) = / d^Jl{x) , (67) 
where the conserved current J§(x) is 

J$(x) = XlgT (r(x)dj(x) - f(x)dj'*(x)) . (68) 
The Hamiltonian for the new field can be constructed in the form 

where the 6-coefficients satisfy the commutation relations 

[b(q),btf)\ = , 
[bHq),tf(q')] = , 

[b(q),bHq')] = 2^^(27r) 3 5(q-q') . (70) 
The Feynman propagator for ip(x) is 

G F (x-x') = t / JL« . (71) 

Assuming that the a-operators and 6-operators commute we can construct a joint system 
and obtain the action 

S = S 4> + S^ > ., (72) 

and the Hamiltonian 

H = H^ + H^ . (73) 

Even though the two theories are individually Lorentz invariant the joint system does not in 
general, as indicated in section 2, exhibit Lorentz invariance. 
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4.1 Mixing Interaction 

In this non-Lorentz invariant context we introduce a mixing interaction between the field 4>(x) 
and ip(x) by adding further term to the action so it becomes 

S = S 4> + S^-w 2 J rf 4 a;(00) a cj)(x)^(x) . (74) 

We have included the factor (S7f2) 2 with the integration weight in order to maintain a formal 
symmetry between the two fields. 

An immediate consequence of including the interaction between the fields is that the theory 
is explicitly non-Lorentz invariant whenever the two lightcones do not coincide. We can still 
carry out the canonical quantization procedure for the joint system. The fields n(x) and co(x) 
conjugate to 4>{x) and ip( x ) respectively are still given by eq(41) and eq(62) and we can still 
impose the standard commutation relations. 

The equations of motion however are now 

n g ^d^d u (f>(x) + n m 2 (j)(x) + w 2 (nn)^(x) = o , 

ng* v dpd v i>(x) + nm 2 4)(x) + w 2 (nn)i<l>(x) = . (75) 
Plane wave solutions are of the form 

W = ( Oft ) = ( M ) = ($) • (76) 



*I>(X) J V fq( X ) J V A l / 

The amplitudes A q and A q satisfy 

n(g^q^q v -m 2 )A q -w 2 {QQ) k -A q = , 

(£'"V Zl , - m 2 ) A q -w 2 {flQ)iA q = . (77) 

In order that these equations be soluble we must impose the determinantal condition 

(g^q^- m 2 )(g^q^-rh 2 )-w 4 = . (78) 

This yields the generalization of the mass-shell conditions appropriate to the model. When the 
mixing term is absent, w 2 = 0, this reduces to the two mass-shell conditions in cq(47) and 
eq(66) . When w 2 7^ and typical components of q^ are large then eq(78) implies that either 
{g^QuQv — m 2 ) or (g^q^qv — m 2 ) is large but not both. If we assume that the second of these 
is large then on writing eq(78) in the form 

4 

id 

(g^q^-m 2 ) = — ^ , (79) 

{g^q^qu - m 2 ) 

we see that there is solution for q^ that approaches the mass-shell of eq(47) with corrections 
that are 0(l/q 2 ) where q is a typical component of q^ . There is an alternative solution in which 
q^ approaches the mass-shell of eq(66) in a similar manner. At large g M we expect that there 
is for each mass-shell a solution that lies near its positive energy branch and another near the 
negative energy branch. That is four solutions in all. By energy in this context we mean p° or 
pP where p a = e^q^ and p a = e^q^ . It clear from the "trivial" case in which the lightcones 
coincide and the theory is Lorentz invariant that there is a restriction on w 2 if all four solutions 
are to have real values of go for all values of qi . That is we must have 

w 4 < to 2 to 2 . (80) 

The marginal case 

W A = m 2 fh 2 , (81) 
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results in a massless mode. Breaking the constraint eq(80) results in solutions with complex 
values for q for sufficiently small values of <fc . The same is true in the Lorentz symmetry 
breaking case. If we impose condition (81) then eq(78) becomes 

(9^q^W T Qaqr)-fh 2 g^q^-m 2 g^q^ = ■ (82) 

This equation clearly has a solution with = . We can expect to find solutions to eq(82) of 
the form q^ = XQ^ that for small A, satisfy 

(mV" + ™ 2 Dw, = • (83) 
That is q^ lies on a light-like mass-shell 

=0 . (84) 

where u — fh 2 / (m 2 +rh 2 ) . The interpolating metric g^ u {u) satisfies g^{Q) = <p" and g^(l) = 
g^ v . It is essential for the causal structure of the theory that g^ u {u) retain a proper lightcone 
structure of the form 

g^{u) = r, ah e\e\ , (85) 

where a is a real vierbein. In particular we require that det g(u) < . Were the determinant 
to vanish and change sign for < u < 1 then no such vierbein could exist and there would be 
a breakdown of causality in the theory at least for values of m and m giving rise to a value of 
u for which det g^(u) > 0. 

As we will see in section 7 this interpolating cotangent metric is of crucial omportance in 
other interacting theories where similar issues of causality breakdown arise. 

4.2 Wavefunctions for the Mixing Model 

If we have two wave function solutions F(x) and F'(x) of the equations of of motion eq(75) with 

F W = ( M ) > (86) 



/(*) 

and similarly for F'(x) then the current J M (x) given by (see eq(48) and eq(68)) 

j»(x) = j;(x) + j;(x) , (87) 

is conserved even though the individual contributions are not. The mixing term in the Lagr- 
nagian leads to a flow of probability between the two parts of the wave function. We can then 
define a scalar product 

(F',F) = J d 3 xJ» , 

that is independent of time. 

If we choose two plane wave solutions as in eq(76) we find for the current 



(88) 



J"{x) - (flA*,A g g^ + tiA q ,A q g^) (q v + ^e"^"^ . (89) 
It is easy to show from eq(75) that, as expected, 

d^(x) = -i(q„-q'JJ»(x) = . (90) 
The scalar product of the two wavefunctions can be obtained from eq(89), thus 

(F q ,,F q ) = (27r) 3 <5(q-q') (ClA*,A q g 0v + Q,A* q , A q g av ) {q v + g^e"^"-^ . (91) 

If (jo and q' are on the same branch of the solutions of eq(78) then q a = q'o when q = q' and 

(F q , , F q ) = (2^) 3 <5(q - q') (CIA* A q g 0v + ClA q A q g°») 2q v . (92) 
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If qo and q' are on different branches of the solutions of eq(78) then for consistency we must 
have 

(W,)=0 • (93) 
This follows from the fact that when q' = q, 

(nA* q ,A q g » + nA*,A q g°»)(q v + ql)=0 , (94) 

which can be proved directly from eq(75) . For definiteness we adopt the following normalization 
convention. Denote the mass-shell based on g^ v as the o-shell that based on g^ v the 6-shcll. 
Then if q lies on the branch that asymptotes the a-shell we use the normalization 

(F q ,,F q ) = 2Clg^q u (2nf6(ci-^) . (95) 

If q lies on the branch that asymptotes the 6-shell then 

(F q ,,F q ) = 2^°V(27r) 3 5(q-q') • (96) 

4.3 Quantization of the Mixing Model 

The standard quantization procedure can be applied to the mixing model. The upshot of course 
is that we can express the quantum fields as a superposition of wave functions with coefficients 
that are annihilation and creation operators thus 

tU ) - / m&FZ«<Wv + / s^^'w + >«■ ■ m 

where the a- and 6-coefficients satisfy the commutation relations in eq(55) and eq(70) and 
commute with one another. The Hamiltonian has the form 

H = 1 wSk^y [q)a{q) + 1 j^Sk^y {q)h{q) ' (98) 

where we have left it implicit that qo takes values on the branch appropriate to the integrand 
in which it appears. 

The Feynman propagators are 

Gfix-x') = -i(O\T(0(x)(f>(x')\Q) 
Gf{x-x') = -i{0\T{<t>{x)i,{x')\0) , 

Gp(x-x') = -i{0\Ty>(x)il>(x?)\0) , (99) 



= Gf(x 


-x' 


Gfix- 


x>) 


Gf{x- 


x') 


ap{x- 


x') 



I f d ' 

Q.J (2^ 



-iq T {x T — x' T ) 



nj (2tt) 4 D{q) 
w 2 f d 4 q e -*iAx T -x' T ) 

my/ 2 J J^y D(qT 

1 f d 4 q (g^q^ - m 2 )e- lc t^ xT - x ' T '> 



nj (2tt) 4 D(q) 



(100) 



where 



D(q) = (g fw q^-m 2 +ie)(g^q^-fh 2 +ie)-w 4 . (101) 

As expected each of the three propagators has in Fourier space, poles in q at points where D(q) 
vanishes. This is of course the determinantal condition in cq(78) . To complete the analysis it is 
necessary to show that these poles lie on the real go-axis and that go-integration contour can be 
threaded through these poles in a way consistent with causality. The discussion of subsection 
4.1 shows that this may not always be possible. In section 5 we consider some special cases that 
illustrate these issues explicitly. 
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5 Special Cases of Bi-Metric Backgrounds 



Coleman and Glashow [8, 17] pointed out that it may be possible to preserve subgroups of 
the Lorentz group even when the full group is no longer an invariance of the theory. These 
possibilities are exemplified by considering the Mixing Interaction model in the following cases. 



5.1 Rotational Invariance 

If we choose the vierbeins so that 

/ Q \ 
10 
e » ~ 10 

V o o o i J 



and 



( tt ^ 

10 

10 

V o o o i / 



(102) 



then 



9nv = 



( n 2 




V o 





-1 









-1 





\ 




-1 J 



and 



( n 2 











\ 





-i 














-1 







V o 








-1 


J 



(103) 



The two lightcones are then 

n 2 (x ) 2 - (x 1 ) 2 - (x 2 ) 2 



3\2 



{X*) 







sf( x y - (x 1 ) 2 - (x z y - (x-y = o 



(104) 



It is obvious that they share invariance under the rotation group in the spatial coordinates. 
The dispersion relation for the mixing model in this case is 



D(q) 



n 2 



2 2 

Qi - ( h 



q\ 



m 



2 2 2 -2 
Ql-Q2- % - m 



W 







(105) 



It is obvious that there are two positive real values for qfi and hence two positive and two 
negative solutions for q provided that w 2 < mm . Explicitly we have 



n 2 n z + n z ji 2 ± J(n 2 ^ 2 - n 2 fi 2 ) 2 + 4n 2 n 2 w* 



(106) 



where fi 2 = q\ + q 2 + q 2 + m 2 and /1 2 
the lower branch, for small q^ has the form 

n 2 n 2 (m 2 + m 2 ) 

n 2 m 2 + £l 2 m 2 



qi + qi+qi+m ■ In the marginal case where w = mm 



{ql + ql + ql) ■ 



(107) 



The implied lightcone structure is perfectly compatible with causality. This is confirmed (rather 
trivially) by computing the determinant of the interpolating metric, we have 



det g^ v {u) = -(utt 2 + (1 - u)Vt 2 ) , 



(108) 



and therefore det g^viu) < for < u < 1 . Finally we note that we can use D(q) in eq(105) in 
eq(100) to construct the propagator with Feynman boundary conditions by setting q 2 — > qfi+ie . 

5.2 Boost Invariance 



If we choose the vierbeins so 


that 











°\ 







a 





e a = 





a 







V 





1 / 



and 



^ 1 ^ 
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a 

\ 1 J 



(109) 
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then f2 = a 2 and £1 = a 2 and 




9 pu> 



/I 





V o 





-a 2 
-a 2 
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-1 / 



and 



9»v 



( 1 
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-a 2 














-a 2 







V o 








-1 


J 



(110) 



The two lightcones are then 



( 2 ; ) 2 -a 2 (x 1 ) 2 -a 2 (x 2 ) 2 -(x 3 ) 2 = . 



(Ill) 



It is obvious that they share invariance under Lorentz boosts in the 3-direction and rotations in 
the (l,2)-plane. Geometrically the two lightcones touch one another in the plane x 1 = x 2 = 
along the lines x° — ±x 3 . 

The dispersion relation for the mixing model in this case is 



1 



2(11+^2) - ll -m 2 



1 



2 (qI + al) - 



m 



w 







(112) 



The solutions for q 2 are 



9o 2 



= 2((^ + ^) ( ^ +<Z22)+2,Z32 + m2 + 



in 




(q 2 + q 2 ) +m 2 — to 2 I + 4w 4 



It is easily checked that detg^u) 



[ua 



(1 - u)a 2 ) 2 < for < u < 1 



(113) 
We therefore 



expect no causality breakdown in this case for any values of m and m 



5.3 Sheared Lightcones 

If we choose the vicrbeins so that 



e , = 



( 1 \ 

10 

10 

\ -K 1 J 



and 



then fl = fl 



9ixu 



( 1 _ 



V 



1 and 
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1 J 
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-K 2 
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— K 








-1 


/ 



The two lightcones are then 



Or ) 2 - (a; 1 ) 2 - (a; 2 ) 2 
(a; ) 2 - (a; 1 ) 2 - (a; 2 ) 2 



9^ 



(x*-KX a f = 
{X* + KX°) 2 = 



(114) 



(115) 



They are obtained from the standard lightcone with n 
positively and negatively respectively. 
The dispersion relation is 

a\2 2 2 2 2\ It \2 2 2 

qo + Kq 3 ) -q 1 -q 2 -q 3 -m ) ((q - nq 3 ) - q 1 - q 2 



(116) 

0, by shearing in the 3-direction 



4 



m 2 ) 



w 



(117) 
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Finding an explicit solution for q of this quartic equation is difficult in general. However in 
the special case for which fh = m the equation reduces to a quadratic equation for with the 
result 

q 2 = (k 2 + l)ql + q\ + q\ + m 2 ± ^ A^q 2 (q 2 + q 2 + q 2 + m 2 ) + . (118) 
If we consider the marginal case w = m we find for small q^ 

dl = Ql + Ql + (1 - « 2 )</3 • (H9) 

For k < 1 (and w 2 < m 2 ) the dispersion relation yields real solutions for q for all values of qi, 
(72 and q% and there arc no problems of causality breakdown. However when k > 1 we obtain 
imaginary solutions for q and we no longer have a proper causal structure for the theory. We 
can confirm this outcome by computing the interpolating metric. We have 



/ 

9p.v = 



-K 2 








(1 - 2u 





-1 














-1 





- 2u)k 








-1 



(120) 



For the equal mass case the relevant value of the interpolating parameter is u = 1/2 . We have 

det g^(u) = -(1 -4u(l -u)k 2 ) . (121) 

Clearly when < k < 1 detg All/ (u) remains negative for < u < 1 . We expect no causality 
breakdown in this case. When k > 1 the determinant will vanish twice in the range < u < 1 . 
When u lies between the two zeros we may encounter causality breakdown. The marginal case 
is k = 1 when a double zero appears at u = 1/2 . If we examine the two lightcones we see that 
for < k < 1 they do share timelike vectors while for n > 1 this is no longer the case. In the 
transitional case with k = 1 the two lightcones touch along a ray that is the x°-axis. In more 
complicated theories than the simple mixing model we will again encounter these issues. 

6 Mixing model with TV Metrics 

We consider briefly a generalisation of the mixing model to the case of N metrics and N fields 
and return to the notation of section 2 . We choose as our mixing interaction a contribution to 
the Lagrangian of the form 

&{x) = \Y^w% J d 4 xV£l^SlW<j> {i \x)<f>( j \x) , (122) 

where Wij — Wji and wu — . The dispersion relation, a generalisation of eq(78), is 

detQ(g) = , (123) 

where the matrix Q(q) is 

Qij(q) = (9 {i)lM/ q^-m 2 i )6 ij -w 2 j . (124) 

Here is the mass parameter associated with field (f>^(x) . The marginal case that allows 
solutions of eq(123) for arbitrarily small values of q is 

detQ(O) = . (125) 

The dispersion relation for small q then becomes 

^m^^ = , (126) 
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where Wi is the i tn diagonal minor of Q(0) . We see that the effective metric controlling the 
small q behaviour of the mass shell is the interpolating metric g^ v where 



with 



5>, 9 W^ , (127) 



" • (128) 



W 1 + W 2 + ■ ■ ■ 

We see again that there must be a constraint on the masses such that g^ v yields a proper 
lightcone. In turn this will impose constraints on the relationships between the individual 
lightcones associated with each of the fields. Of course a detailed analysis of the general case 
is rather complex. Nevertheless the above discussion shows the importance of the interpolating 
metric in understanding the causal structure of the theory. The interesting question is to what 
extent these issues re-emerge in a study of more general models of interacting fields. In the next 
section we will examine such a model for the case of two quantum fields. 



7 Interacting Field Theory - Quartic Interaction 

We can construct a model that includes include scattering effects by allowing quartic interactions 
for the fields. A model that is relatively simple but non-trivial has an action of the form 



where 



with 



and 



with 



and 



S = S 4 , + Si, + S I , (129) 
5 = f (FxtlC^x) , (130) 



C^x)= l -{g^d^{x)d v ^x)-rnlcj> 2 {x)) , (131) 



J d n xnC^{x) . (132) 



£j>( x ) = \ ^{x)d v <f>{x) - mlip 2 (x)) , (133) 



with 



Si = J d n xd{x) , (134) 

£ / ( a; ) = -^fi (^)) 4 -^n (^(x)) 4 -^(fi n )5(^)) 2 (^)) 2 . (135) 

In the above the squared bare masses for the fields cf>(x) and ip(x) are, m 2 , and to 2 , respectively. 
We have also introduced the bare quartic coupling constants A , A and a together with bare 
metrics gfo u and g^ v since we will find that these latter are also renormalised. The bare vierbein 
determinants fi and Oo ar e related to the bare metrics in an obvious way. Our aim here is to 
investigate the renormalisation of the bare parameters in low order of perturbation theory using 
a dimensional regularization scheme [18]. We have therefore expressed the action of the theory 
in terms of n-dimensional integrals. 
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7.1 Feynman Rules 

We first formulate the perturbation series interms of the bare parameters, expanding in powers 
of A , Ao amd a ■ The Feynman Rules for computing the terms in the series for the Green's 
functions are essentially the same as for the conventional theory (in which the metrics coincide) 
except that the vertices and propagators incorporate extra factors involving powers of l^o and 

n ■ 

In constructing Feynman diagrams we associate the propagation of a 0-field with a solid line 
and a factor (see section 4) 



iG(q) 



1 



g^l^lv -ml + ie 



where is the momentum flowing through the line. Similarly the propagation of a ^>-field we 
associate with a dashed line and a factor 

*G(q) = ' 



The vertices in a Feynman diagram associated with the 4> 4 , ip 4 and (jptp 2 interactions are 
— iAofio, — *AoOo and — iao(fio^o) 1 ^ 2 respectively. Of course we enforce momentum conservation 
at each vertex and finally integrate over loop momentum with a weight 

d n k 



(27t)" 

Standard symmetry factors are applied as appropriate to each graph. 
7.2 Four Point 1PI- Amplitudes 

To examine the renormalisation properties of the sreies we compute the lPI-amplitudes derived 
from the truncated Green's functions. We will carry out the analysis to second order in the 
coupling constants. 

The lowest order contributions to thelPI-amplitude for four external </>-lines, M4, are in- 
dicated in Fig 2. Since we are interested only in the divergences of these contributions it is 
sufficient to evaluate them with zero momentum on the external lines. Taking into account 
the symmetry factors for the loop diagrams and the fact that there are three distinct ways for 
attaching the external lines in each case we find 

l M 4 = -zA ^ + ^(^A ) 2 /+^H ( 7 ) 2 ^J , (136) 
Z Z S Zq 



where 



and 



r d n k 

~ J (2^F 



[ V _ -^0 r f n \ (n-4) nr7 x 



dnk ( 1 \ _ T ( 2 n )m {n - 4) (138) 

(27r)» I, g^kph -ml + ie) ~ (4tt)"/ 2 \ 2 ) ' [ M) 



Confining attention to the divergences at n = 4 this yields 
iM4 — — zJ7q 



Ao+ (4^ (A ° + ff ° 2) ^ 1 4 



(139) 



The corresponding results for the lPI-amplitude for four external zero momentum ^-lines is 
M4 where 

_ [_ 3 _ 1 

iM 4 = -iCl A + — -r^CXl + al) 

(47T y n — 



(140) 
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Figure 2: Feynman Diagrams for the amplitude M4 



The lowest order contributions to the lPI-amplitude with two external 0-lines and two external 
^-lines, M 2i 2, are indicated in Fig 3. Again we set the external lines to zero momentim we have 

m-i\o)(-i<Jo)(no/n ) 1 / 2 j + 2(-i<j ) 2 K , [ ' 

where I and J arc defined above and K is given by 



l k i 

r) n g^kftK -m% + ie g^ v k^k v -m% + ie 

We have not given the evaluation of / and J explicitly because it is entirely the same as the 
corresponding integrals in standard theory. However the evaluation of K introduces a new 
feature involving the interpolating metric g^ v (x) = xg^ v + (1 — x)gQ V . We therefore do provide 
an explicit evaluation. Using the standard representations of the Feynman propagators 



9Q V k ll k v -ml + ie 
i 



ie) 



poo 

= / dX ie iX ^(9o U k^K-ml + 
poo 

5o kpkv ~ m o + %e Jo 



(143) 



Making the change of integration variables Ai, A 2 — > x, A where Ai = Ax and A 2 = A(l — x) we 
find 

pi poo p mi, 

K= dx dXXe -iX(xm 2 +(l-xm 2 -ie) / _±_^_ ( ,iXg^(x)k l ,k 1/ _ 

Jo Jo J (27r) n 

In order to evaluate the /c-integral it is necessary to require that g^"{x) does not become singular 
for < x < 1 . As we have explained in section 3.3 this requirement is guaranteed if we demand 
that the two lightconcs do overlap in a manner that permits the existence of coordinate rays 
that are timelike in both metrics. We therfore assume this to be true from now on. We then 
obtain 



K=-±- tdx [ X d \\e-^ m2 o+^-^ 2 o-^—J^=(?-) n/2 
(2tt)" J Jo ^/-detg^(x) VzA7 



that is 



(4tt)"/ 2 Jo y/-detg^ v (x) 

We are concerned only with the divergence so we may simplify this to 



(145) 



(146) 



dx 

(47r) 2 n-4i ^J- dct.g^(x) 

We have then 

T 2 



^ 2 , a = -i(nono)V»L + ^-i- + 7 fL^-) , (us) 

\ (47r) 2 n — 4 (4-7r) 2 n — 4/ 
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X X>- X'-< >°< 

Figure 3: Feynman Diagrams for the amplitude M22 

where 



dx 

We can now determine the renormalisation structure for the coupling constants. We expand 
the bare parameters to second order in rcnormalised parameters by setting 



,4— n 



,4—n 



a 


n 




4 




a 




n 




4 




b 




n 




4 



where 



a = ^- n [°+z—A X +~ X ) + —^ 2 ) > ( 15 °) 



3 

(4tt)2 ' 
4 



6 = - 



(47T) 



2 



c = -(Sy ■ (151) 

We assume, as will be clear later, that the renormalisation structure for the bare metrics has the 
form g^" = g^ v + 0(\ 2 , a 2 ) and g^" — g^ v + 0(A 2 , a 2 ) so in the above calculation we are free to 
replace all metric quantities by their renormaliscd versions. It then follows that the amplitudes 
M4, M4 and M 2j 2 are finite to second order in the renormalised parameters. 

A significant result to emerge from the above analysis is that the the lightcone structure 
influences the details of the renormalisation procedure. However if we were to allow the two 
metrics to coincide then the above expansion for the bare couplings in terms of the renormalized 
couplings would coicidc with the standard result for the theoretical model we are investigating. 
If we allow the two metrics to approach the situation in which their timelike overlap reduces to 
zero we would find that the parameter 7 diverges therefore the renormalisation procedure would 
exhibit a singularity in the expansion of the bare parameter ctq interms of the renormaliscd 
parameters. The precise physical significance of this is unclear but we obviate the problem by 
postulating that the two metrics must retain a timelike overlap for the quantum field theory to 
exhibit a proper causal structure. This is in line with our conclusions from investigating the 
simple mixing model in section 4.1. 

7.3 Two Point lPI-Amplitudes 

We denote the sum over lPI-amplitudes with two external 0-lines by E(g) . This amplitude 
contributes to the inverse two point Green's function thus 

G2 1 (q) = {l (g^q ft q„-ml) + X(q) . (152) 

Similarly for the lPI-amplitudes with two external V'-hnes we have 

G^(q) = n Q (g^q^-fh 2 Q ) + t(q) . (153) 
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Figure 4: One loop bubble Feynman Diagrams for the amplitude 



7.3.1 One Loop Bubble Diagrams 

The Feynman diagrams up to one loop, appropriate to computing iH(q) are shown in Fig 4. We 
have 

■w v = -^Ao f d n k i -i* /Q\ 1/2 f d n k i 

l ^ q> 2 J (2nrgZ»q^-mt+ie + 2 {ClJ J (2tt)™ - m§ + ie ' ^ 

That is 

= -iOo^^ra-n^Xm 2 )"/ 2 - 1 -^^^^ . (155) 

A parallel calculation yields for the lPI-amplitude for two external ^-lines 

iZ(q) = -ffi ^^r(l-n/2)(m 2 )"/ 2 - 1 -i(O O ) 1 / 2 ^^r . (156) 

Working to first order in the renormalised masses and coupling constants and retaining only 
the divergent pole terms at n = 4 we find 

ts( g ) = -in 2 1 . m 2 - i(nn)va ^ i ^2 _ (157) 

(47r) 2 n - 4 (47r) 2 n - 4 

Here we have introduced the renormalised mass parameters m 2 and m 2 which are related to the 
bare masses by the expansion to first order 

ml = m 2 (l + -A_A) + ^ 2 (-^. 

,2 f, , d t\ , 2 Z' / 



ntf = m^l + —Xj +m ^-^aj . (158) 

We have also replaced Oo and J7o by their renormalised values n and fl since the difference 
involves terms of second order in the couplings. If we choose 

d = d = 1 



(4vr) 2 

/ = 



1 /f^ 1/2 



(4tt) 2 \n 

1 1/2 



> = "wUJ ' (159) 

then we see from eq(152) and eq(153) that the poles at n — 4 in the expressions for GJ 1 ^) and 
G ! 2~ 1 (</)cancel as required by the dimensional regularization procedure. 

For future reference we now use eq(150) and eq(151) to expand the right sides of eq(155) 
and eq(156) to second order in couplings. Dropping non-pole terms we list the coefficients of 
the powers of couplng constants for iT,(q) as follows 

m. 1 

A : -in 



(47r) 2 n - 4 
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a : -i(fifi) 



-in 
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1 1 
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6 1 1 

(4tt) 2 V(^-4) 2 + n-42 
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47T/J : 
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1/2 _ 
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1 1 



(4tt) 2 \(n-4) 2 n-42 
Here ^(z) is the logarithmic derivative of T(z) . 



(6 + J) (log 



777' 
47T/J ; 



-V>(i)-i) + d 



(160) 



7.3.2 Two Loop Bubble Diagrams 

The two loop bubble diagrams that contribute iS(g) are shown in Fig 5. Because they are 
already second order in the renormalised couplings we can replace all bare quantities by their 
lowest order expansion in renormalised parameters. The loop integrals are easily evaluated and 
yield the following results where we list the coefficients of the powers of coupling constants 
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(161) 



If we drop finite contributions and retain only the poles at 77 = 4 we obtain 
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(162) 
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Figure 5: Two loop bubble Feynman Diagrams for the amplitude £((/) 





Figure 6: Two loop Rising Sun Feynman Diagrams for the amplitude £((/) 



7.3.3 Two Loop Sunrise Diagrams 

The two loop sunrise diagrams contributing to iS(g) are shown in Fig 6. Their evaluation is 
performed in appendix B. The hrst diagram is 0(A 2 ) and gives rise to the pole terms 
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(163) 



This is the same result for general linear coordinates as in the standard calculation for 4 -theory. 
The second diagram a similar but more complex outcome. It is convenient to present it in three 
separate pieces. The first piece which depends on m 2 , is 



iE(g) 



Here the quantity £ is 
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1 (4tt) 4 
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The second piece which depends on to 2 is 
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o y/-dctg» v {u) Jo 1 - « V v /-det.g^(u) 
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(164) 
(165) 

(166) 
(167) 
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The third term yields the 0{q )-dependent term, 
where 

The matrix d^ v is, from its derivation, symmetric. 
7.3.4 Renormalisation of the Metrics 

The contribution to iT,(q) in eq(169) together with the term 0(q 2 ) in cq(163) induce 
malisation of the bare metric g^ v . We set 

This implies that 

n = n(i-l x2 ^ fTP + ^ hTP ) . 



2 n-4 J 

Now to 0(q 2 ) we can write, ignoring non-pole contributions, 

G- 2 \q) = G;\0) + Qotfw " ^^^ frfr - ^n^/^ " 
Using cq(170) and eq(171) we can remove the poles from the 0(q 2 ) term by requiring 

7 ~~ 12(47r) 4i/ ' 

and 

It is convenient to define x = 9nv 

6^ so that 

1 f 1 duu „ 

x= ^y v-det^H^^ ■ 

Following a parallel calculation we have for the other metric 



and 



We find 



and 



where 



V 2 / 



J 12(4^ 



^/-det <?^0) 
It is convenient to define x = g M „ 6^ so that 

1 f 1 du(l - u) 



~2tt/ V-det^H^^ 9 " 
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7.3.5 Mass Renormalisation 



In order to complete the mass renormalisation to second order in the coupling constants we first 
add up all the appropriate contributions to iT,(0) . Referring back to eq(160), cq(162), eq(164), 
eq(166), and eq(168) we find that the combined result for the coefficients of the powers of the 
renormalised coupling is 
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We then require that to 2 , is a series in the renormalised coupling constants with coefficients that 
are poles in n at n = 4 in such a way that f^o^o ~ ^(Oj ni 2 ,fh 2 ) remains finite at n = 4 . To 
second order in the couplings this yields 



ml = m 2 M + m 2 



1/2 



M 



(183) 



where 
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A 2 



(4tt) 2 n-A (4tt) 4 V(«-4) 2 
a 2 / 2 

+ 



+ 



1 1 



3n-4 



(l-C- X /4) \ 
(4tt) 4 V(«-4) 2 n-4 J 



and 



M 



1 (7 

+ 
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(184) 



(185) 



The important point is that there is no residual dependence on [i in the expressions for 
the bare masses. It is clear however that as was the case with the coupling constants, the 
renormalisation procedure involves the metrics, through the interpolating metric as well as 
directly. 



7.4 Renormalisation Group 

Strictly speaking we have checked only that the renormalisation procedure works to second order 
in the couplings. Although this is a non-trivial result we should, in principle, exhibit the result 
to all orders. This is, as with all field theories, a formidable problem made more difficult in our 
model by the manner in which the metrics enter the process. We shall nevertheless assume the 
correctness of the renormalisation procedure and use it to deduce the the renormalisation group 
equations for the parameters of the theory. 
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The renormalisation group equations are essentially the statement that the bare parame- 
ters of the theory are invariant under changes of the scale parameter \i . We have then the 
requirement 

d , d — d ,„o^ 

(i— A = ^ — A = ^i — cr = . (186) 
dfi d\x d\i 

In addition in our case we must require 

<%# = "°- (187) 

If we introduce the variable t — log(/i///i) where /ii is a reference scale then from eq(150) and 
eq(151) we find to second order in the coupling constants, the results 

A = (n-4)A +(4 ^(A 2 + ( x 2 ) , 

X = (n-4)A+^(A 2 + ( x 2 ) , 

a = (n-4)a+^a(A + A) + ^a 2 , (188) 

where i-derivatives are indicated by an over dot. From eq(170), eq(173), eq(174), eq(176), 
cq(178) and cq(179) we find again to second order 

»" " <£r?' + • < 189) 

Eqs(188) and eqs(189) are inter-related because 7 depends on the metrics. 

In order to investigate the structure of these equations it is useful to set g^ v — Aj^ and 
g^ v = Bj^ u where det 7^ = det ^ u = — 1 . The equations for the evolution of A and B become 

6(4tt) 4 2(4tt) 4 

A 2 a 2 x 
6(4tt) 4 + 2(4tt) 4 



B I,:/)!, ' "l'; )/>' , (190) 



The equations for 7^ and 7^" are 



cr 2 



V"' = (yV^ - 46^) 

7 2(4^) 4lX7 J 



f = 2^ ( ^"" 4 ^ ) ' (191) 

where W = 6^ v /A and 6^ = ff^/B . It turns out that 7, x, X, W and depend only 
on 7 M!/ and . From its definition in cq(149) we can show using the change of integration 
variable 

*'= *A + l-*)B • (192) 

that 

7= / = , (193) 
7o V- dot 7"" (a;') 

where 

f"V) = xY" + (1 - ^Ot"" ■ (194) 
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Similarly we have 

1 f 1 dx'x' 

ijf ^CPgr ^'^ ■ (195) 

where 7 M ^(a;') is the inverse of 7'"' (a;') . We also have 



@^ _ 1 f 1 dx'(l-x') 



X = 

X = • (196) 

We can reduce the renormalisation group equations further by expressing them in a frame 
corresponding to the vierbein of one of the metrics. We introduce a and E^ a such that 

V" = E» a E\rf b , 

f v = ^ a E\jf h . (197) 



The two vierbcins arc related thus 



E\ = E\N\ , 



E\ = E\N b a , (198) 

where N b a is the inverse of N b a . It follows that 

N\ = E\E\ . (199) 

Since 7 M1/ is invariant under Lorentz transformations of the vierbein, that is E^ a — > E^ b L b a , 
we cannot immediately obtain an equation for E^ a unless we impose a further constraint. The 
relevant constraint in this case is to require that 

W a r? b E\ = E» a r? h W h . (200) 

The point is of course that were the constraint not to be satisfied it is always possible to construct 
a Lorentz transformation L a b (t) that will ensure that eq(200) is satisfied by the transformed 
vierbein. Similar remarks apply to the construction of E^ a . It follows that 

r v = 2E\ V ab E\ , 

*T = 2E\ V ab E\ , (201) 
It then follows that eq(191) implies that 

E\ = ^y 4 ( X E\-^E% Vac ) , 

E\ = ^^{xE\-^E% Vac ) , (202) 
We can compute the renormalisation group equation for N a b from cq(202) and the result 

N b a = E»^-E\N c a + E\) . (203) 

We have then 

N b a = ^4 ((X -x)N b a + 4&^E b ^E d v r, dc N c a AQ^E b ^ vVda ) . (204) 
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It is convenient to introduce the matrix fj ab (u) and its inverse r) a b{u) where 

f] ab (u) =ur 1 ab + (l-u)r 1 cd N a c N b d , (205) 

so that 

r"(u)=E\E\r b (u) . (206) 
It follows that det7 M!/ (u) = det?7 ab (w) . The parameter 7 can be expressed as 



du 

y/-dct fj ab (u) 



Similarly we find 



2 Jo y/- det f] ab {u) 
1 f du(l-u) , 



* = 9 / / y d flcd{u) , (208) 

2 7o v/-dctr/ ab (ii) 

We find also 

z Jo v det V \ u ) 

e^E b d E\ Vda = \ f f U{ \~^ ( , n bq f, qP {u)N\ . (209) 
A Jo V — det r] ab (u) 

From eqs(204) and (209) we can show that 

4 l 47r J Jo v — det rj ab [u) 

where 

M 9<, (u) = (1 - u)ri qd - uN" r N* tf" . (211) 

7.5 Fixed Points of the Renormalisation Group Equations 

It is clear from the results of subsection 7.4 that the parameters A, A, a and N a b satisfy a closed 
set of differential equations. The evolution of the field renormalisation variables A and B is 
determined subsequently. An important aspect of the renormalisation group equations is the 
structure of the fixed points. 

We can easily check that the right side of eq(210) vanishes when N b a — 6 b . More generally 
this is true when N b a is a Lorentz transformation. This is to be expected since the theory then 
has only the standard lightcone structure and is Lorentz invariant. The neighbourhood of the 
this fixed point in iV-space can be understood by setting N b a = S b + n b a where we regard n b a 
as infinitesimal. From eq((210) we find to first order in n b a 

n b a = lj^(n b a + n c dV bd Vac ) . (212) 

It follows that the difference of the two terms on the right is time independent. This corresponds 
to the fact that eq(210) is invariant under Lorentz transformations of N b a from the right. We 
will simplify the situation and assume that n b a lies in a subspace for which the two terms are 
equal. We have then 

n b a = ~n b a - (213) 
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The solution is 

n\(i)=n\(0)expj^VgpJ . (214) 

If then a were to approach a non-zero fixed point value as t — > — oo we could conclude the 
theory would acquire Lorentz invariance in the infrared limit. In fact we will see below that the 
situation is somewhat more complicated. 

The analysis of the fixed point structure in the near Lorentz invariant limit is simplified by 
noting that to first order in n b a the parameter 7 = 1 and the renormalisation group equations, 
eq(188), for the coupling constants reduce to the standard form for the Lorentz invariant model. 
If we adopt the Wilson approach [19] and work in n = 4 — e dimensions then we see that there 
are a number of fixed points in coupling constant space in the limit of weak Lorentz symmetry 
breaking. 

The most obvious fixed point, A, lies at the origin of coupling constant space. It is clear 
from cq(188) that for e > this fixed point is infrared unstable and correspondingly ultraviolet 
stable. There is a second fixed point, B, with a = and 

A = A=^e = 2r, , (215) 

where 7% = ((47r) 2 /6)e. We can show that this point is infrared stable and ultraviolet unstable 
with respect to all three coupling constants by examining the renormalisation group equations 
in its neighbourhood. Set A = 2r st + u and A = 2r* + w and treat u, w and a as small. From 
cq(188) we find 

6r* 



(4tt) 2 

W = (4^ 
10r* 



u (4tt) 2 " ' 

learly u, w and a all vanish in the infrared limit. In particular 



(216) 



a=aiexp {^4 ' (2i7) 

where <j\ is the value of a at t = . 

There is a third fixed point, C, in the domain of positive coupling constants that lies in the 
plane A = A and satisfies the equations 

A 2 + a- 2 - 2r*A = , 

\ X +\°~ r * =° • ( 218 ) 

The solution in the appropriate sector of positive couplings is 

A = A = A, = 0.36117% 

a = <7* = 0.77187% (219) 

By examining the neighbourhood of fixed point C we can easily check that it is ultraviolet 
stable in one direction and infrared stable in two other directions. The situation can be largely 
understood from Fig 7 which exhibits the three fixed points lying in the plane A = A . The 
interesting renormalisation group flows lie along the sides and in the interior of the "triangle" 
ABC . The flow along side BC is ultrviolet stable at C the flow along the side CA is infrared 
stable. The other flows originate in B and and end in A . The arrow on each flow indicates the 
direction of increasing t . If we set the masses m and fa to zero for simplicity, we can think of 
each flow as representing a theory. 
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Figure 7: Renormalisation Group flows of the coupling costants, for A = A. The fixed points are 
A (ultraviolet stable), B (infrared stable) and C (ultraviolet stable in the SC-direction, infrared 
stable in the CA direction). 

A BA-flow corresponds to a theory for which the ultraviolet behaviour is controlled by A 
and the infrared behaviour by B . We can therefore assert that such a theory is aymptotically 
free and at large scales will behave as if the effective couplings are A = A = 2r* and u = . 
That is at large distances the two helds <j>(x) and if>(x) will interact only weakly. From eq(214) 
we have 

„» o ( t )=„» o (0)exp{ i g^-«p{^-l}} . (220) 
In the infrared limit we see that the vanishing of a at the fixed point implies that 

n\(0=n\(0)cxp{- i5 J )2r J . (221) 

It follows that because a vanishes in the infrared limit as indicated in eq(217) the two lightcones 
remain distinct. However to the extent that the two fields (p(x) and ip(x) are decoupled and 
each field is invariant under its own version of the Lorentz group. Similar comments apply in 
the ultraviolet limit. At intermediate values of t the breakdown of Lorentz invariance shows up 
becoming maximal when a acquires its largesr value on the trajectory. 

The theory corresponding to the CT-flow is ultraviolet free but in the infrared limit is 
controlled by the couplings A* and <r* corresponding to the fixed point C . Such a theory will 
exhibit a breakdown of Lorentz invariance at large distances. The theory corresponding to the 
BC flow will be controlled in the ultraviolet limit by the parameters of the fixed point C and 
hence will exhibit Lorentz breakdown in this limit. In the infrared limit the theory will again 
be controlled by the parameter values at B as discussed above. 

If we pick a theory corresponding to an interior JSA-flow that starts out from B close to 
the BC-flow then it will remain close to the BC-fiow for a large range of t and will spend a 
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long "time" in the neighbourhood of C thus mimicking the BC-theory to high energies before 
departing from C and moving close to the CA-flow. The corresponding theory therefore appears 
to be controlled at high energies by the fixed point C but at very high energies reconstitutes 
itself as a CA-theory. The intervening stage that is dominated by C exhibits maximal Lorentz 
breakdown. Of course the B A- flow along the lower edge of the "triangle" ABC corresponds to 
(7 = and hence no interaction between the two fields. 

When these mass parameters are not zero the discussion of the infrared limit will require 
some modification to take this into account. We do not deal with this in this paper. Furthermore 
when we return to four dimensions, that is e — > 0, eq(221) implies that finally n b a vanishes and 
the two lightcones coincide. In this limit also the fixed point values of A and A vanish and the 
theory becomes trivial as expected for a quartic scalar theory in four dimensions, as well as 
Lorentz invariant. However in the spirit of the e-expansion [19] we might hope that our analysis 
applies at least qualitatively to three-dimensional scalar theories. It would be of great interest 
to check this directly by other means. 

It is worth noting that the infrared stable fixed point C does not require 7=1 for its 
existence. The discussion of the infrared limit in this case therefore holds good even when the 
separation of the lightcones is finite. 



7.6 Finite Departures from Lorentz Symmetry 

To investigate the possibility of other nontrivial fixed points in general is not easy. We present 
two simple cases which can be explored relatively easily. 



7.7 Rotational Invariance 

In the first of these we assume that the matrix 7V a b relating the two vierbeins frames in theory 
has the form 
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(222) 



This represents a situation where the theory remains invariant under rotations about the 3-axis 
and under boosts along the 3-axis. Eq(210) implies 



OJ 



a 



1 + uj 4 oJ 4 log(a; 4 ) 

l-OJ* (1-w 4 ) 2 



(223) 



It is easily checked that the coefficient of a 1 only vanishes when oj = 1 . This is therefore only 
fixed point value for oj . If we set oj = 1 + p and assume p is small then we find 



3 (4tt) 



(224) 



which is consistent with eq(213) . In this sector of the theory therefore we find only the fixed 
point structure already encountered. 



7.8 Sheared Lightcones 

It is of interest also to exmine a situation in which the relation between the two metrics of the 
theory incorporates a shear of one lightcone relative to the other. In order to achieve this in a 
way that is consistent with the renormalisation group equation (213) we set 



N a , 



c 

1 

1 

-as 



-s/a \ 



c J 



(225) 
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where c = coshi/> and s = sinhtp for some hyperbolic angle ip . The interpolating metric fj ab (u) 
is given by 



( u+(l-u)(c 2 -s 2 /a 2 ) -(l-u)(a-l/a)sc \ 

-10 

0-1 

\ -(l-u)(a-l/a)sc u + (1 - u)(c 2 - a 2 s 2 ) J 



(226) 



There is no real loss of generality if we assume that ip is positive and a > 1 . We have 
det7y ab = - A where 

A = 1 - u(l - u)(a - l/a) 2 s 2 . (227) 

This is the form of the determinant for the interpolating metric that arises when the shear 
k = (a — l/a)s . It is obvious that we must impose the constraint n < 2 in order to avoid the 
interpolating matrix fj ab (u) becoming singular for < u < 1 . The integrand on the right side 
contains the matrix eq(213) 



r] pq (u)M qd (u) 



(A B \ 

(l-2u) 
(l-2u) 



V c 











D ) 



where 
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B 
C 
D 



= ^[-u 2 (c 2 - s 2 /a 2 ) + (1 - u) 2 (c 2 - a 2 s 2 )] 
= ^(u 2 + (l-u) 2 )(a-l/a) S c , 
= -^(u 2 + (l-u) 2 )(a-l/a) S c , 
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A 



[-u 2 (c 2 - a 2 s 2 ) + (1 - u) 2 (c 2 - S 2 /a 2 )] 



(228) 

(229) 
(230) 
(231) 
(232) 



Using the symmetry A under the interchange ao(l-ii)we find 

du 



o \P~ det i) ab {u) 



% q {u)M^ d {u) 
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(233) 



where 



with 



P = -R(a 2 - l/a 2 )s 2 
Q = 2R(a-l/a)sc , 



R 



-f 

Jo 



duu 



A 3 / 2 ' 

The renormalisation group equations for a and ip can be inferred from eq(213). They are 

T 2 



(234) 
(235) 
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ri?(a-l/a) 2 (l + s 2 ) . 



(236) 



These equations imply a renormalisation group equation for the shear parameter k namely 

T 2 



K = 
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(M 



Rk(4- k 2 ) 



(237) 
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It is obvious from eq(235) that R depends only on k so that eq(237) involves only the coupling 
constant a and n . Superficially there appears to be fixed point for n at k = 2 . However it is 
clear that R has a divergence for this value of k and the leading singularities can be computed 

as 

as k — V 2 from below. The presence of the pole at k = 2 removes the zero in the right side of 
cq(235) and hence any possibility of a fixed point. Eq(235) then implies that in the infrared 
limit the shear parameter will renormalise down to zero and the theory will again approach 
the fixed point in which the lightcones coincide. The presence of the logarithmic singularity at 
k = 2 means that for k > 2 the renormalization group equations of the theory cease to be real. 
A similar divergence is exhibited by the parameter 7 we have 

(239) 

This underlines the point, already evident in the calculation of the perturbation series for 
the Green's functions that the the two lightcones must be overlap in the sense of there being 
spacetime directions that are timelike in both metrics. 

The indications from these two examples are that there are no infrared stable fixed points 
other than those already discussed. To complete the proof of this it would be necessary to 
extend the analysis to the case where N a b takes on arbitrary values in the group SL(A) though 
this is a considerable task and beyond the range of this elementary approach. 



8 Conclusions 

We have investigated simple models of quantum field theory in a multimetric background. Our 
starting point was to seek a causal structure for such theories by requiring that the various 
metrics in the model shared a foliation of spacetime that was spacelike with respect to each of 
the metrics simultaneously. This suggested the relevance of a parametrisation of the relation- 
ship between metrics that involved a rescaling of the coordinates together with an appropriate 
shearing transformation that preserved the spacelike foliation. 

In the case of a bi-metric theory with metrics g^ v and g^ v and two scalar fields without 
interaction our initial restriction on the metrics was sufficient to permit the combining of the 
two quantum fields into a single quantum system. Formally invariance under Lorentz trans- 
formations is lost although the evolution each field is uninfluenced by this. When a simple 
transition interaction is introduced the breakdown of Lorentz invariance is real since each field 
becomes aware of the other. Although formally the theory can still be formulated as a causal 
combined quantum system in fact there are circumstances that can bring about a breakdown of 
causality. In assessing the presence of acausality we are led to construct the interpolating metric, 
g^(u) = ugV + (1 — u)g^ Vl with < u < 1 . For certain choices of interaction strength this 
metric, for particular values of u, controls the small momentum behaviour of the dispersion re- 
lation. Even though the individual metrics give rise to conventional light cones the interpolating 
metric may become singular for u in its standard range. There will under these circumstances 
be a range of u for which the interpolating metric corresponds to a pseudo-Lorentz metric with 
diagonal entries (+1, +1, —1, —1) . This is incompatible with a causally propagating system. In 
order to avoid this possibility it is necessary to further restrict the relationship of the two metrics 
so that the interpolating metric is never singular for u in its standard range. We show that the 
relevant restriction is that the lightcones of the two metrics overlap and that therefore there 
are directions in spacetime that are timelike in both metrics. When the model is generalised 
to encompass more than two metrics the same analysis reveals the relevance of an interpolating 
metric that is a positive superposition of all the metrics in the model. We have not investi- 
gated this situation in detail but the general point remains that to avoid potential acausality 
we should require the configuration of metrics to be such that the interpolating metric remains 
non-singular. 
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A bi-metric model in which the fields have quartic interactions is studied. We analyse 
the perturbation series for the model to second order in the coupling constants. Dimensional 
rcgularisation is used to deal with the renormalisation of the theory. Two issues stand out. 
First, the interpolating metric g^ u {u), again plays a crucial role. The evaluation of the Feynman 
diagrams is unambiguous only if g^{u) remains non-singular throughout the standard range 
for u . It follows that we again must impose the requirement of overlapping lightcones in 
order to have a causal theory. Second, given this restriction, it is possible to carry through 
the renormalisation program, at least to second order. However this involves not only the 
standard renormalisations of couplings, masses and fields but also the two metrics themselves. 
We therefore must begin by formulating the theory with appropriate bare metrics and obtain, 
scale dependent, renormalised metrics as part of the calculation. 

The renormalisation procedure gives rise to renormalisation group equations for the renor- 
malised quantities. The equations for the renormalised couplings and masses are similar to those 
that would obtain were the metrics to coincide. They are however influenced by the relationship 
of the two metrics. In turn the renormalisation group equations for the metrics are of course 
dependent on the renormalised couplings. We can separate off the dependence of the overall 
scale of the metrics and identify it with what is conventionally viewed as field renormalisations. 
We then show that we are free to view the theory from the vierbein basis of the first met- 
ric which then takes the form rf h while the second metric takes the form r\ — N a c N b d rf d . 
The renormalisation group equations then condense a little and involve only the couplings the 
masses, and N a b . The matrix N a b involves the scaling and shearing transformations discussed 
in section 3 . 

When N a b is near the identity the renormalisation group equations for the couplings reduce 
to those for the Lorentz invariant case. Following Wilson and Kogut [19] we work in 4 — e 
dimensions and are able to pick out significant fixed points of the these equations that allow us 
to identify theories that are well controlled at scales in the ultraviolet and infrared limits. Of 
course the well known triviality of quartic scalar theories re-emerges when e — > and all fixed 
points return to the origin. However in the spirit of the e-expansion [19] we might argue that 
the set of theories is an approximation to a corresponding set of theories in three dimensions 
obtained by setting e = 1 . We examine two situations in which N a b departs strongly from the 
identity. One involves scaling and the other a shearing. In neither case does it seem likely that 
there are further fixed points with large Lorentz symmetry breakdown. The analysis is however 
not complete. 

It is worth remarking that if we generalise the quartic model in an appropriate way to involve 
more metrics and fields we would expect to reproduce the same pattern of restrictions in the 
lightcones of the metrics by requiring any two light cones to overlap and share in the manner 
explained a set of timelike directions. This is a pairwise condition and does not mean that all 
metrics share a timelike directions. However the result emerges from second order perturbation 
theory calculations and may be modified at higher order. Intuition suggests that it is likely 
that we will find ultimately that all participating metrics must share in a mutual overlap of 
spacetime vectors that are timelike in each of them. In these circumstances we would expect 
that notwithstanding the breakdown of Lorentz invariance there are time variables that permit 
the introduction of a Wick rotation and a Euclidean formulation of the theory. 

It would be interesting to test these ideas in higher order and in more physical theories 
such as the Lorentz non- invariant extensions of QED and the Standard Model [9, 10, 11, 20] 
where the interest lies in the energy scale for the onset of Lorentz symmetry breakdown and the 
behaviour of the theory above and below this scale. 
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A Singularities of Interpolating Metrics 



Although not completely general it is sufficient for the purposes of demonstrating the result to 
make the choice for the two metrics as folows 



=M\M\r, aT , 



where 



M", 



( 1 


M 

V o 

where M is a general real 3x3 matrix. We have then 
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-M T M 
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(242) 



It is convenient to choose a coordinate basis that diagonalises M 7 M. leaving g^ v in the form 
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( 1 \ 
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\ -1/c 2 J 



(243) 



The associated lightcone is 
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The other metric is chosen to be 
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wherc S a ,, is a shear transformation of the form 
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(244) 
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The intersections of the two lightcones with the plane x a = 1 are respectively the ellipsoid 
centred at the origin 



+ 



+ 



= 1 



and the sphere centred at the point (a; 1 , x 2 , x 3 ) = (vi,v 2 ,v 3 ) 

(x - Vl ) 2 + (x 2 - v 2 f + (x 3 - v 3 ) 2 = 1 . 



(249) 



(250) 
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When the centre of the sphere is sufficiently near the origin the sphere and the ellipsoid will 
overlap. The points in the interior of the overlap correspond to coordinate rays that are timelike 
with respect to both metrics. When displacement of the sphere is sufficiently large there will be 
no overlap with the ellipsoid and the two metrics will not share any timelike rays. The marginal 
case which is of interest here is one in which the sphere lies outside the ellipsoid but remains 
touching it. At the point of contact the coordinates satisfy both eq(249) and eq(250). There is 
a further condition that the normals to the surfaces are anti-parallel, that is 



(251) 



for some A > . We can use eq(251) to eliminate the coordinates from eq(249) and eq(250) 
with the results 
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An appropriate combination of eq(252) and eq(253) yields 
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These equations determine A and impose a constraint on (^1,^2,^3) 
is g^„(u) = ug^ + (1 - u)g^ v , that is 



( 1 - uw 2 

uv 2 
uv 3 



It follows that 

where 

F = 1 - u(l 



-(u + (l-u)/a 2 ) 





1 



uv 2 


-(u+(l-u)/b 2 ) 




The interpolating metric 
\ 



uv 3 



-(«+(!- 



det g^ = -[u + 



u + 



b 2 



u + 



a 2 u + (l-u)/a 2 b 2 u+(l-u)/b 2 c 2 u + (l-u)/c 2 



The interpolating metric becomes singular when F vanishes. That is when 



1 



1 



1 



a 2 I + {I - u) / {ua 2 ) b 2 1 
If we make the identification 



(l-u)/(ub 2 ) c 2 l + (l-u)/{uc 2 ) 



(255) 



(256) 



(257) 



(258) 



(259) 



then eq(258) becomes identical to eq(254). The marginal situation that we are investigating is 
that F has coincident zeros in u . We ensure this by demanding in addition that 



dF 
du 







(260) 



This immediately yields eq(253) after imposing cq(259). Finally we note that eq(259) implies 
that the range A > corresponds to < u < 1 and therefore the requirement of antiparallcl 
normals corresponds to the appropriate range for u . 
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B Rising Sun Diagrams 



We exhibit the calculation for both of the Rising Sun diagrams in Fig 6. Because these diagrams 
are of second order in the coupling constants we can replace bare quantities by their lowest order 
renormalised expressions. 



B.0.1 Three Internal </>-lines 

It is useful to discuss the calculation of the first diagram in a certain amount of detail. It provides 
guidance in how to deal with the second less conventional diagram. Using the Feynman rules 
we find that the first diagram yields a contribution to iT,(q) given by 

*£(«) = -^AV) 4_n ^(«V) , (261) 

where 

2 [°° d n k l( l n k 2 i i i 

(q,m ' J (2tt) 2 » D(k 1 ,m 2 )D(k 2 ,m 2 )D(k 1 + k 2 -q,m 2 ) ' [ ' 

with 

D{k,m 2 ) = g^k^k v -m 2 + ie . (263) 
From the Schwinger representation of the propagator, 



J dX ie iXlD< ~ kl ' m2) etc , (264) 



D(k\, m?) 
we have 

/d n k\ d n k 2 f 
(2 ^2n J dXid\ 2 d\ 3 

enp{i[X 1 D{k 1 ,m 2 ) + X 2 D{k 2 ,m 2 ) + X 3 D(k 1 +k 2 -q,m 2 )}} . (265) 

We introduce the change of variables Ai = Xxi, X 2 = Xx 2 and A3 = Xx 3 where x\+x 2 +x 3 = 1 . 
We have then 

L(q,m 2 ) = J J d Xl dx 2 dx 3 5(l - Xl -x 2 - x 3 ) J dXX 2 e^ E , (266) 

where 

E = Xig^k^k^ + x 2 g^ v k 2 ^k 2v + x 3 g^(k! + k 2 - g) AI (fc 1 + k 2 - q) v - m 2 + ie . (267) 
We displace the origin of the integration variables to the stationary point of E by setting 

; X2X 3 , 7^ 

ki = -. , , ^q + Ki , 

(XlX 2 +x 2 x 3 +X 3 Xl) 

k 2 = -, ,q + K 2 . (268) 

(XiX 2 + X 2 X 3 + X 3 Xi) 

The exponent E takes the form 

E = £ + - ;™ s g^q,q v - m 2 +ze , (269) 

{XxX 2 + X 2 X 3 + X 3 Xl) 

where 

£ = x l9 ^K^K lv + x 2 g^K 2 ^K 2v + x 3 g^{K x + + K 2 ) v . (270) 

The momentum integrations can now be performed 

ggigga <Ag = _ Jl 1 1 

(2tt) 2 ™ (4Tr) n (iX) n ( Xl x 2 + x 2 x 3 + x 3 X!) n / 2 ' k ' 
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Since our aim is to calculate the ultraviolet poles at n — 4 wc can expand in powers of q to 
second order. We obtain, incorporating the above evaluations, 

L(q,m 2 ) = L(0,m 2 )+g^q^M{m 2 ) , (272) 

where 

Tin 2n • ^ tvq m 2yn- 3 f dx 1 dx 2 dx 3 S{l - X! - x 2 - x 3 ) 

L(0,m ) = —i-, — —rl(3 — n)(m ) / ; : — j- — , 273) 

1 ' (4tt) 4 1 A ' J (x lX2 + x 2 x 3 + x 3Xl ) n / 2 K 1 

and 

2 = _ 2 4 f dx^dx^l-x, x 2 x 3 )x,x 2 x 3 

(47TJ 4 J (#1#2 + X 2 X 3 + X 3 Xl) n ' A+L 

As explained by Ramond, the evaluation as it stands, of the integral in eq(273) is tricky. We 
achieve a more tractable case by first introducing into the integrand the trivial factor (x\ +x 2 + 
x 3 ), which evaluates to 1 against the 5-function. We then use the symmetry between xi, x 2 and 
x 3 to replace this with 3x 3 . This procedure corresponds to the use of an integration-by-parts 
identity introduced by 't Hooft and Veltman. We have then 



Q 2 

L(0, m") = — z- 

wherc 



B _ r dxidx 2 dx 3 S(l -xi-x 2 - x 3 )x 3 
J (xix 2 +x 2 x 3 +x 3 xi) n / 2 

To complete the evaluation we eliminate the <5-function by performing the a^-integration, then 
make the change of variables X\ = xy and x 2 = (1 — x)y (the range of all variables is (0, 1)). 
We have 



where 

' \ZV i 

(1 - y(l - x(l - x))) n / 2 



B' = J dxdyy 1 -"/ 2 

For our purposes it is sufficient to evaluate B' at n = 4 . We find B' = 1 and hence 



(278) 



B= ^+3 , (279) 

n — 4 



and as a result 



L(0,m 2 ) = 6i ^ (fl _ 4)a (l-(n-4)(#) + 3/2)) . (280) 
From eq(274) we see that M(m 2 ) has only a simple pole at n = 4 . The residue is contained in 

where 



dxidx 2 dx 3 S(l — xi — x 2 — x 3 )xiX 2 x 3 



(282) 



Using the change of variables explained above we find 

A = [ dxdy - X{1 r X){ }- y \^ = I • (283) 

J U (1 - y(l - z(l - z))) 3 2 v ; 

Hence 

y ' 2(4tt) 4 n-4 V ; 
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Finally combining these results we find that the pole contributions to iT,(q) are contained in the 
result 



i£(g) = -in 



A 2 m 2 / m 2 



n—4 r 



l-(n-4) U(l) + 



(4tt) 4 V4tt^ 2 
This gives immediately the result in eq(163) 
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B.0.2 Two Internal ^-lines and One Internal </>-line 

Applying the Feynman rules to the second diagram in Fig 6 yields a contribution to iT,(q) given 

by 

iE( ? ) = -^V) 4 - n iL(g,m 2 ,m 2 ) - ( 2 86) 

where 

r/ 2 -2\_ f d n kid n k 2 i i i . . 

L[q,m,m)-J { ^ )2n ^ ^ ^ + ^ _ ~ > (287) 

where 

D{k,rh 2 ) = g tlv k li k v -m 2 + ie . (288) 

Introducing the Schwinger representation for the propagators and making the appropriate 
changes of integration variables along the same lines as the previous calculation we obtain 

L(q, m 2 ,m 2 ) = J / da^da^l - an - x 2 - x 3 ) J d\\ 2 e lXE , (289) 

where 

E = xxg^ 'k llM k 2v + x 2 g» l 'k 2)X k 2v 

+x 3 g^(k 1 +k 2 -q)^(k 1 + k 2 -q) v -(x 1 +x 3 )fh 2 -x 2 m 2 + ie . (290) 

The stationary point of E is given by the requirements 

dE _ dE 
dk lfl dk 2fl 

That is 



. (291) 



xig' lv k lv + x 3 g' lv {k 1 + k 2 -q) v = , 

x^k^ + xsg^ih + h-q), = . (292) 

These equations have the solution k\ = kf^ and k 2 = k^ where 

= X2X3 g, u (u)g»°(u)q a , (293) 

p x\x 2 + x 2 x 3 + x 3 xi 

and 

*4? = x ^ g^? T 9 r P {u)g p »q v . (294) 

Note the appearance of the interpolating metric g^{u) . The argument u is 

u = . (295) 

XlX 2 + x 2 x 3 + X 3 Xl 

Also 

l_ u= ^ , (296) 

X\X 2 + X 2 X 3 + X 3 X 1 

We can then write 

C=(l-^fSrpMA • (297) 
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Using the identity 

{ugr + {\-u)gr)g av {u)=% > ( 2 98) 

we can show that 

k ( $=q li -ug lip g' n 'q v • (299) 

We then have 

(k^+k^-q),= g w rq u . (300) 

XlX 2 + X 2 X 3 + X 3 Xl 

We shift the origin of the momentum integration variables by setting 

k2 tl = k { V+K 2ll . (301) 

We have then 

E = £ +£ R - (xi + x 3 )m 2 - x 2 m 2 + ie , (302) 

where 

£ = xigrK ilx Ki v + x 2 g^K 2 ^K 2v + x 3 g^{Ki + K 2 )„{K X + K 2 ) v , (303) 
and the remainder, £r, is of course the value of E at the stationary point and is given by 

£r = x^gTkf^ + x 2 r v kf^ + x 3 g^(k{ 0) + fcf - g)„(fc<°> + fcf - q) v . (304) 

Note that £r = 0(q 2 ) for small q . 

After performing the momentum integrations we have 



/ 



WV»m)-.. ! __J = . (305) 

(2tt) 2 ™ (4tt)™ v (x!X 2 + x 2 x 3 + x 3 X!)"/ 2 y/- det g^(u) 



Since we only wish to compute the poles at n = 4 in L(q, m 2 7 m 2 ) it is sufficient to consider the 
low q expansion 

L(q,m 2 ,m 2 ) = L(0,m 2 ,fh 2 ) + M^q^q u . (306) 
We obtain these terms by expanding 

e i\e R = 1 + iX£r ; (307) 

which is correct to 0(q 2 ) . Performing the A integration we find 

2 _ 2 . Q f dx\dx 2 dx 3 8{l — x\ — x 2 — x 3 ) 



2 _ 2 . n f dxidx 2 dx 3 8(l - xi - x 

L(0,m ,m ) = — T(3-n) / . =- 

V 7 (4tt)« v ; i V-det^H 

[(xi + x 3 )m 2 + x 2 m 2 ] n ~ 3 



and 



(xix 2 + X 2 X 3 + X 3 X!)™/ 2 

g?x 2 g?x 3 <5(1 — xi — x 2 — x 3 ) 



(308) 



^/- det (^"(u) 

[(xi + x 3 )m 2 + i 2 m 2 ]"" 4 , 
(xix 2 + x 2 x 3 + x 3 xi)™/ 2 



(309) 



The evaluation of £r is as follows. Using the results for k^ and fc 2 °^ we find 
£r = x ™ [u-g^g^(u)g Ta q a g^q p 

XiX 2 +X 2 X 3 +X 3 Xi 



(1 - u)g» T g Te g™q a g„ a g aX g x0 (u)g^q p ] . (310) 
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Using the identity in eq(298) this becomes 

X1X2X3 



-R 



X\X2 + X2X3 + X3X1 



{9 aT grM~9 W )Qaq P ■ (311) 



The matrix (jf g T /j(«)<? Alp ) can be shown explicitly to be symmetric. 

It is convenient to re-express L(0, to 2 , to 2 ) as a sum of two terms in the form 

L(0,TO 2 ,TO 2 ) = -i^ i r(3-n)[m 2 J L 1 +TO 2 L 2 ] , (312) 



where 



( dxidx2dx 3 S(l — x\ — x 2 — Xs) 



= i 



and 



yj- det g^{u) 

(xi + x 3 )[(xi + x 3 )m 2 + x 2 to 2 ]"" 4 
(xxx 2 + x 2 x 3 + X3X1)"/ 2 

dx\dx2dx 3 8{l — xi — x 2 — x 3 ) 
y/- det g» v (u) 

X2 [{x 1 + x 3 )fh 2 + X2m 2 ] n ~ 4 

(X1X2 + X 2 X 3 + X3X1)"/ 2 



(313) 



(314) 



We now eliminate the (5-functions by integrating over x 2 and introduce the change of variables 
x\ = xy and X3 = (1 — x)y with the result 

(315) 



1 — y(l — x(l — x) 
and 

= r dxdyy 2 ~ n / 2 [ym 2 + (1 - ?/)m 2 ]"~ 4 
1 7 V- det ¥^(1 - y(l - z(l - a;))) n/2 

We invert eq(315) to replace y by u as the integration variable and obtain, ignoring contributions 
that vanish when n = 4, 

Li - [m) J T^mm ' ( } 



with the result 



_ y Mi-*?"™ . (318) 

r(4-n) v ; 7 ^-dct.^^) v ; 



By a parallel calculation we find 



(r(n/2-l)) 2 n _ 4 [duu(l-uf-l 2 



L 2 = ^(to 2 )™- 4 / ' =r- . (319) 

r(n-2) V ; 7 ^/-det^H 

We have also 

Combining these results with eq(286) and retaining only the pole terms at n — 4 we obtain the 
results in eq(164), eq(166) and eq(168). 
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